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ABSTRACT 


A  classical  result  of  Gilbarg  states  that  a  simple  shock  wave  solution  of 
Euler '8  equations  is  compressive  if  and  only  if  a  corresponding  shock  layer 
solution  of  the  Navier-Stokes  equations  exists,  assuming,  among  other  things, 
that  the  equation  of  state  is  convex.  An  '"entropy  condition*  appropriate  for 

T' 

weeding  out  ""unphysical  shocks  in  the  nonconvex  case  has  been  introduced  by 
Ts— FT* Liu.  For  shocks  satisfying  his  entropy  condition,  Liu  showed  that 
purely  viscous  shock  layers  exist  (with  zero  heat  conduction).  Dropping  the 

convexity  assumption,  but  retaining  many  other  reasonable  restrictions  on  the 

f'lZQo 

equation  of  state,  «e  ^construction  example  of  a  (large  amplitude)  shock  which 
satisfies  Liu's ' entropy  condition  but  for  which  a  shock  layer  does  not  exist 

Pq. 

if  heat  conduction  dominates  viscosity. ^We  also  gives^a  simple  restriction, 
weaker  them  convexity,  which  does  guarantee  that  shocks  which  satisfy  Liu's 
entropy  condition  always  admit  shock  layers.  — — — *■ 
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SIGNIFICANCE  AND  EXPLANATION 


A  traditional  method  of  determining  which  shocks  (which  are  discontinuous 
solutions  in  gas  dynamics  or  elasticity)  are  "physical"  is  to  ask  whether  the 
discontinuity  can  be  obtained  in  the  limit  of  vanishing  viscosity  and  heat 
conduction*  For  a  constant  profile  shock,  one  then  seeks  corresponding 
constant  profile  "shock  layer"  solutions  when  viscosity  and  heat  conduction 
are  nonzero.  If  the  equation  of  state  is  convex,  a  simple  condition 
distinguishes  those  discontinuities  which  admit  shock  layers:  The  shock  must 
be  compressive,  or  entropy  should  increase  across  the  jump.  Convexity  does 
not  hold  for  some  materials,  however.  In  the  nonconvex  case,  Liu  proposed  a 
revised  "entropy  condition"  to  select  physical  shocks.  For  small  amplitude 
shocks,  the  present  author  has  shown  that  a  shock  layer  always  exists  if  Liu's 
entropy  condition  is  satisfied.  Here  we  present  an  example  to  show  that  this 
is  not  the  case  for  strong  shocks,  unless  the  equation  of  state  satisfies  some 
extra  restrictions.  We  also  supply  some  extra  restrictions  which  are  shown  to 
be  appropriate. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


NONEXISTENCE  OP  A  SHOCK  LAYER  IN  GAS  DYNAMICS 
WITH  A  NONCONVEX  EQUATION  OP  STATE 

Robert  L.  Pego 


f 1.  Introduction 

In  1951,  D.  Gilbarg  [2]  showed  that  for  each  stationary  shock  wave 
solution  of  Euler's  equations  of  gas  dynamics  in  one  space  dimension,  there 
exists  a  corresponding  smooth  solution  (called  a  shock  layer  or  shock  profile) 
of  the  Navier-Stokes  equations  with  viscosity  and  heat  conduction,  provided 
that  the  thermodynamic  equation  of  state  satisfies  a  short  list  of 
restrictions  given  by  Weyl  [11],  including  a  convexity  condition*  The  shock 
wave  is  required  to  satisfy  the  Rankine-Hugoniot  jump  relations,  and  the 
entropy  condition,  which  ensures  that  entropy  increases  along  particle  paths 
in  Euler's  equations.  Under  Weyl's  restrictions,  the  entropy  condition  simply 
requires  shocks  to  be  compressive,  i.e.,  density  must  increase  along  particle 
paths . 

More  recently,  interest  has  developed  in  equations  of  state  which  may  be 
nonconvex  ( [5] , [9] ,  ( 10] ) .  T.-P.  Liu  introduced  an  entropy  condition 

appropriate  for  shocks  in  this  situation.  He  showed  that  a  shock  satisfies 
his  entropy  condition  if  and  only  if  a  purely  viscous  shock  layer  exists, 
i.e.,  a  shock  layer  for  the  equations  with  zero  heat  conduction  [4]. 

The  purpose  of  this  note  is  to  present  an  example  concerning  restrictions 
on  the  equation  of  state  under  which  Liu's  entropy  condition  does  select  just 
those  shocks  for  which  shock  layers  exist  for  any  positive  values  of  viscosity 
and  heat  conduction.  We  explicitly  construct  an  equation  of  state  such  that 
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for  hm  shook  wavs  which  satisfies  Liu's  entropy  condition*  no  corresponding 
shook  layer  can  exist  if  the  heat  conduction  dominates  the  viscosity.  (A 
shock  layer  does  exist  if  heat  conduction  is  small  compared  to  viscosity.) 

The  equation  of  state  satisfies  all  of  Wsyl's  restrictions  except  the 
convexity  condition*  and  satisfies  other  conditions  which  have  appeared  in 
recent  work  on  the  Riemann  problem  in  gas  dynamics  ([5] ,[10]).  He  can  give  a 
simple  reetriction  weaker  than  convexity  (satisfied  by  van  dar  Meals  gases  in 
regions  of  hyperbolicity,  for  example)  under  which  all  shocks  satisfying  Liu's 
entropy  condition  always  admit  corresponding  shock  layers.  (see  Theorem  1.) 

An  example  such  as  we  describe  may  be  of  interest  in  one  dimensional 
thermo (visco) elasticity,  where  for  quite  general  constitutive  relations  it  is 
desirable  to  determine  the  smoothing  effect  of  viscous  and  thermal  dissipation 
(e.g.  see  [1]).  Nonconvex  equations  of  state  also  occur  for  materials 
exhibiting  phase  transitions.  In  this  paper,  however,  we  require  that  Euler's 
equations  remain  hyperbolic  in  the  region  of  interest.  Effects  other  than 
viscosity  and  heat  conduction  must  be  considered  to  determine  "structure"  in 
phase  transition  zones;  see  [8]  and  [9] . 

The  Navier-Stokes  equations  in  one  space  dimension,  written  in  Lagrangian 
form,  are 

Tt‘  \  -  ° 

d.D  ut+  Ph  ’  (wVT)h 

&t+(pu)h  -  <»«VT)h  +  <xVT)h 

Here  h  denotes  the  Lagrangian  mass  coordinate,  t  is  time,  t  is  specific 
volume,  u  is  velocity,  p  is  pressure,  8  is  temperature,  &  is  energy  desity  per 
unit  mass,  and  u  and  X  are  the  coefficients  of  viscosity  and  heat  conduction, 
respectively.  4  ■  e  +  u2/2  ,  where  e  is  the  internal  energy  density  per  unit 


mass.  We  assume  that  t  and  8  determine  the  thermodynamic  state  of  the 
material,  and  that  e,  p,  and  the  entropy  S  are  given  by  sufficiently  smooth 
equations  of  state  e  *  e(T,0),  p  ■  p(t,0),  S  *  S(t,0)  .  These  functions  are 
related  through  the  Gibbs  relation 

(1.2)  0  dS  -  de  +  p  dt 


The  coefficients  u  and  X  may  also  depend  smoothly  on  t  and  0.  The  quantities 
U,  X,  t,  6,  and  p  are  positive.  For  convenience,  we  shall  denote  the  triple 
(t,  u,  & )  by  U,  and  the  thermodynamic  state  by  Z,  specified  by  the 
pair  (t,0)  (or  by  (T,p),  see  assumption  (1.10)). 

Euler's  equations  are  obtained  from  (1.1)  by  setting  u  »  X  •  0  .  A  shock 
wave  traveling  at  speed  s  is  a  weak  solution  of  these  equations  of  the  form 


(1.3) 


U(h,t) 


U_  for  h  <  st 
U+  for  h  >  st 


The  two  end  states  must  satisfy  the  Rankine-Hugoniot  relations 


0 

*  -S(T  -  T  ) 

-  (u  - 

u  _) 

+  - 

+ 

(1.4) 

0 

=  -s(u+  -  u_) 

+ 

+ 

1 

pj 

0 

-  -s<a+  -  &_) 

+  <p+u+ 

-  P_u_) 

We  call  any  solution  (1.3)  satisfying  (1.4)  a  "simple  jump  solution",  or 
"jump".  It  is  customary  to  reserve  the  term  "shock"  for  jumps  which  satisfy 
the  entropy  condition. 

What  determines  a  jump?  The  speed  s  may  be  determined  up  to  sign  from 
any  suitable  pair  of  thermodynamic  states  Z+  and  Z_  by  the  relation 
(1.5)  0  -  s2(t+  -  x_)  +  <p+  -  p_) 

In  our  analysis  it  suffices  to  consider  only  back-facing  shocks,  for  which 
s  <  0.  (Particle  paths  in  (1.1)  are  vertical.)  Now  u+  -  u_  is 
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determined.  Imposing  the  third  equation*  we  recall:  The  states  Z+  and  Z_ 
determine  a  simple  jutqp  solution  with  s  <  0  if  and  only  if  the  Hugoniot 
function  vanishes: 


H(Z+,Z_) 


s+  -  e_  +  V^p+  +  p_)(x+  *  T_)  * 


Fixing  ZQ,  the  set  of  states  Z  satisfying  H(Z*Z^>  ■  0  is  typically  a  curve* 
called  the  Hugoniot  curve  with  center  Z^  in  the  (x,p)  plane. 

A  shock  layer  for  the  junqp  (1.3)  is  a  traveling  wave  solution  U(h-st)  of 
(1.1)  with  0(C)  ♦  U  as  £  ♦  Plugging  in  and  integrating  once,  U(£) 

must  satisfy 

0  -  -S(T  -  X  )  -  (U  -  M±) 

uu^/x  -  -s(u  -  u±)  +  (p  -  p^) 

wuu^/t  +  AQ^/x  -  -s(4  -  4±>  +  pu  -  P±“t 

The  first  equation  determines  u(£)  from  x(£),  so  we  obtain 


(1.6) 


AS  /T  -  -8(8,  -  8  +  (T  -  X  )  (p  +  SU>>  5  L(  X, 0) 

s  X  XX 


UX^/X 


-  x±>  +  P  -  p±> 


=  M(x,0) 


These  equations  yield  an  autonomous  system  of  ODEs  in  the  (x*0)  plane.  It  is 
clear  that  Z+  and  Z_  are  rest  points  for  this  system.  A  shock  layer  exists 
for  the  jump  (1.2)  if  and  only  if  the  system  (1.5)  admits  a  trajectory  Z(£) 
connecting  Z_  to  Z+,  i.e.,  Z(C)  ♦  Z  as  £  ♦ 

He  remark  that  equations  (1.6)  are  identical  to  those  used  by  Gilbarg  [2] 
for  the  stationary  shock  layer  in  Eulerian  form,  if  we  introduce  an  eulerian 
space  variable  x  *  /^x(n)dn  and  recognize  that  -s  is  the  mass  flux 
pu  «  const  ,  where  o  -  1/x  is  the  density. 


In  order  to  motivate  Liu's  entropy  condition  for  the  jump  (1.3)  and 


examine  the  implications  of  nonconvexity*  we  briefly  consider  the  simple 
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example  of  isotherm  l  gas  dynamics  in  Lagrangian  form*  The  equations  are 

\  •  “h  •  0 

"t  tp<T)h  '  <wh/T>h 

A  simple  jump  with  speed  s  connecting  (x_,u_)  *nd  <T+»U+)  ^or  M  ®  satisfies 

0  -  s2(x+  -  t_)  +  p(r+)  -  p(T_) 

A  shock  profile  <x,u)<£),  with  £  -  h-st,  must  satisfy  (eliminating  u) 

(1.7)  V T^/T  *  _^(S2(T  -  T±)  +  p(t)  -  p(t±))  =  M(T> 

Assuming  that  p(T)  is  decreasing,  but  possibly  nonconvex,  any  two  states 
t  ,  t  determine  a  jump  with  speed  s  <  0*  A  shock  layer  exists  for  that  jump 
if  and  only  if:  sgn  M(x)  “  sgn  (x+  -  x_)  for  t  between  x+  and  x_  . 

This  is  the  "entropy  condition"  which  selects  those  jumps  which  admit  a  shock 
layer  for  these  equations.  The  situation  is  pictured  in  Fig.  1,  where  the 
phase  portrait  of  (1.7)  is  indicated  on  the  t  axis  of  the  graph  of  p(x). 
Observe  that  since  p  is  not  convex,  rarefaction  shocks  can  exist  admitting 
shock  layers,  e.g.  x  -  x_  <  X.  ■  X,  .  Also,  a  compressive  jump  need  not 
admit  a  shock  layer,  e.g.  x_  ■  x  >  x+  *  T1  . 


Figure  1.  Shock  profiles  in  isothermal  gas  dynamics 


For  the  nonisothermal  Euler's  equations,  we  now  state  Liu's  strict 
entropy  condition  £or  a  jump  (1.3)  determined  by  Z+/  z_  with  s  <  0  : 
the  condition  is  that 

sgn  M(Z)  «  sgn(T+  -  T_)  for  all  Z  *  (T,p) 
s(E)  between  Z+  and  Z_  on  the  Hugoniot  curve 

with  center  Z+  ( or  with  center  Z_ ) . 


With  s  <  0,  this  condition  simply  means  that  the  line  M(r,p)  * 

s2(t  -  r  )  +  p  -  p  ■  0  should  lie  above  (below)  the  Hugoniot  curve  between 
z  z 

Z+  and  z_  if  t  >  t+  (t_  <  t+).  The  restrictions  we  shall  impose  on  the 
equation  of  state  imply  that  with  some  choice  of  center,  the  Hugoniot  curve  is 
regular  between  Z+  and  Z_,  so  the  condition  s(E)  is  well  defined. 

The  restrictions  we  place  on  the  equation  of  state  are  as  follows: 

(1.8) 

This  means  that  heat  conduction  is  a  dissipative  effect  in  (1.1)  for  X  >  0. 
(1.9) 

Euler's  equations  are  strictly  hyperbolic  just  when  this  holds. 


eQ (t ,0 )  >  0 


p  (T,S)  =  dp/dt  <  0 

T  S  const 


(1.10) 


Pq(t,0)  >  0 


Thus  the  map  (t,8)  +  (T,p(T,0))  is  one  to  one  on  the  domain  of  states  £)  . 
The  image  of  this  map  we  designate  We  assume 


0  is  convex. 
P 


(1.11) 

These  restrictions  above  are  the  same  as  those  imposed  by  Gilbarg  [2] , 

following  Weyl  (11],  except  that  we  omit  the  convexity  assumption 

d2p/dt2  >  0  . 

S  const 

Under  these  restrictions,  the  Hugoniot  curve  is  described  by  the 
following.  The  proof  is  deferred  to  §2. 

Proposition  1 .  Fix  Zg  *  ^Tg'Pg)  *-n  ^p*  Th®h  in  the  region  t  >  tg,  the  set 
{  Z  e  Qp|  H(Z,Zg)  *  0  }  consists  solely  of  a  monotonicaily  decreasing  curve 
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p  *  hg(t)  passing  through  Zq. 

We  now  describe  our  positive  results  concerning  existence  and  uniqueness 
of  the  shock  layer  for  any  given  u,  X,  positive  functions  of  t  and  6. 

Theorem  1.  Assume  ( 1 .8)-( 1 • 1 1 ) •  Also  assume 
<1.12)  eT<T,8)  >0  in  Q0 

Suppose  Z+,  Z_  in  satisfy  H(Z+,Z_)  ■  0  and  (1.13)  below.  Then  the  simple 
jump  solution  determined  by  Z+,  Z_  and  s  <  0  from  (1.5)  admits  a  unique  shock 
layer  if  and  only  if  the  strict  entropy  condition  s(E)  is  satisfied. 

In  the  absence  of  (1.12),  compressive  shocks  ( t_  >  t+  if  s  <  0)  always 
admit  a  shock  layer: 

Theorem  2.  Assume  (1.8)-(1.11).  Suppose  Z+,  z_  in  satisfy  H(Z+,Z_)  -  0 
and  (1.13)  below,  and  suppose  t  >  t+.  Then  the  conclusion  of  Theorem  1 
holds. 

Weak  shocks  ( |U+  -  U_|  small)  always  admit  a  shock  layer,  under  the  most 
basic  of  restrictions: 

Theorem  3.  Assume  only  (1.8), (1.9).  Pix  ZQ  in  G0.  Then  if  Z+  and  Z_  are 
sufficiently  close  to  ZQ  and  satisfy  H(Z+>Z_)  »  0  ,  the  simple  junp  solution 
determined  by  Z+,  Z_  with  s  <  0  admits  a  shock  layer  (unique  in  a  neighborhood 
of  Z0)  if  and  only  if  the  strict  entropy  condition  s(E)  is  satisfied. 

Theorem  3  is  proved  in  [7]  as  a  consequence  of  a  general  theorem  on 
singular  viscosity  matrices  for  systems  of  conservation  laws.  The  proof  of 
Theorems  1  and  2  is  extremely  similar  to  that  of  Giibarg  [2] ,  and  will  be 
presented  briefly  in  $2.  The  extra  condition  imposed  is  a  technical  "domain" 
condition, 

(1.13)  The  line  segment  lies  in  Q0  , 

nee<* >d  so  certe 'n  constructions  are  valid.  (A  similar  assumption  is  made 
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tacitly  by  Gilbarg. ) 

Our  main  result  is  a  "couterexaraple" : 

Theorem  4.  There  exists  an  equation  of  state,  defined  in  a  domain 
satisfying  ( 1 .8)-( 1 . 11 ) ,  and  a  pair  of  states  Z+,  Z_  which  satisfy 

H(Z+,Z_)  »  0  and  (1.13),  such  that  the  jump  solution  determined  by  Z+,  Z_ 
and  s  <  0  satisfies  the  strict  entropy  condition  s(E),  but 

-  if  X/p  is  sufficiently  large,  no  corresponding  shock  layer  exists, 

-  if  X/y  is  sufficiently  small,  a  unique  shock  layer  exists. 

From  Theorem  2,  this  shock  must  be  a  rarefaction  shock,  t_  <  T+  .  We 
stress,  however,  that  Liu's  result  shows  that  a  purely  viscous  shock  layer 
(with  X  -  0)  exists,  so  entropy  does  increase  across  this  shock,  i.e., 

S(Z_)  <  S(z+)  . 

This  equation  of  state  must  violate  (1.12)  and  the  convexity  condition. 
However,  the  following  conditions  can  be  satisfied: 

(1.14)  eT(x,p>  >  0 

(1.15)  Pt(t,9)  <  0 

Condition  (1.14)  (along  with  ( 1 .8)-( 1 . 10 ) ,  basically)  was  imposed  by  Liu  [5] 
to  guarantee  that  the  Riemann  problem  for  gas  dynamics  has  a  unique  solution 
globally.  Condition  (1.15),  along  with  (1.8),  implies  (1.9)  and  implies  that 
constant  states  for  the  Navier-Stokes  equations  are  linearly  stable  (see  [6], 
[7]). 

Finally,  we  briefly  remark  that  although  we  have  presented  this  example 
in  the  context  of  gas  dynamics,  it  is  valid  in  the  context  of  ID  thermovisco¬ 
elasticity  as  well.  One  may  redefine  the  equation  of  state  outside  a  compact 
domain  so  that  the  restrictions  imposed  by  Dafermos  [1]  are  satisfied.  The 
key  step  is  to  appropriately  specify  the  function  p^(T,6)  globally  (e.g.  see 


section  3) 


§2.  Existence  and  uniqueness  of  the  shock  layer 


We  begin  this  section  by  gathering  several  useful  thermodynamic 
calculations,  then  we  proceed  to  prove  Theorems  1  and  2.  Take  t  and  0  as 
independent  variaoles.  Using  (1.2),  and  equating  mixed  partials  of  S(t,0), 
we  obtain  the  standard  identities 

(2.1)  0  SQ (i,8)  -  eQ(T,9) 

(2.2)  0  S  (i,0)  -  e  (r,0)  +  p  -  0po(T,0) 

Considering  t  and  p  independent,  we  get 

(2.3)  8  S  (t,p)  »  e  ( t , p )  *  e„(T,0)  0  (i,p) 

p  p  op 

(2.4)  0  ST(t,p)  -  eT(t,p)  +  p 

and  also 

(2.5)  p  <T,S)  =  -S  (T,p)  /  S  (T,p) 

T  T  P 

The  assumptions  (1.8)-(1.10)  imply  that  the  derivatives 

So(x,0),  s  (t,6),  S  (t,p)  and  S  (i,p)  are  all  positive. 

0  t  p  t 

The  following  identity  helps  relate  the  entropy  condition  s(E)  to  the 
ODEs  (1.6). 

Proposition  2.  Fix  Z+  and  Z_  with  H(Z+,Z_)  *  0.  Let  s  <  0  satisfy  (1.5). 
Now  L(t,8)  and  M(i,0)  in  (1.6)  are  determined,  and  for  any  Z  *  (1,9)  we  have 

(2.6)  H(Z,Z  )  *  -L(Z)  +  V^T  ~  T  )(-sM)(Z) 

z  ~s  ± 

It  follows  that 

H(Z,Z+)  -  H(Z,Z_)  -  V^1t_  -  T+)(-sM)(Z) 

Proof.  H(Z,Z^)  +  L(Z)/s  *  -  u^)  +  (i  -  T±)(V^p  +  p^)  “  p^  +  su) 

*  (T  -  T^)(1/y(u  +  U±)  +V^p  -  p±)  -  su)  «  Vj(T  -  1^)(-SM)(Z) 
(Recall  that  u  -  u±  =■  -s(i  -  i^).  ) 

Proof  of  Proposition  1.  With  ZQ  fixed,  H(Z,ZQ)  is  a  function  of  Z  -  (i,p). 
Using  (2.2),  Hr<i,p)  -  0St<i,p)  +  V^PQ  -  p)  >  0  if  p  <  pQ,  and 
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H  <T,p)  -  •  (T,p)  +V^T  -  T  )  >  0  if  T  >  t  . 
p  p  *  U  0 

It  is  convenient  now  to  calculate 

(2.7)  _~Ii0  ( 1 , 8 )  -  e0(T,0)  _^Lt(t,8)  -  8p0(t,0)  +  »M(t,0) 

(2.8)  -sM0(t,8)  -  p» ( x » 8 )  -sM  (x,0)  -  a2  +  p  (x,  0 ) 

0  0  T  T 

To  get  ,  use  Proposition  2  and  (2.2). 


Proof  of  Theorem  2.  Let  Z+,  Z_  in  fl0  satisfy  H(Z+,Z_)  *  0,  (1.13),  and 
t_  >  t+/  and  let  s  <  0  satisfy  (1.5).  Suppose  first  that  the  entropy 
condition  s(E)  is  satisfied.  Our  proof,  like  Gilbarg's,  is  based  on  an 
analysis  of  the  curves  (L  ■  0)  and  (M  -  0)  which  pass  through  Z+  and  Z_. 
Since  L0  and  MQ  are  positive,  0  is  a  given  function  of  t  along  each  of  these 
curves,  0  »  Hr),  8  «■  m(T>  respectively.  Vte  plain  that  t(x)  and  m(x)  are 


defined  for  x+  <  x  <  x  and  satisfy 
i)  Hr)  <  m(T) 


ii)  r<t)  <  0 


for  T  <  T  <  T 
for  X+  <  T  <  X_ 


Indeed,  the  set  (M  ■  0)  in  is  a  straight  line  segment,  by  (1.11).  Pulling 

back  to  the  (t,0)  plane,  m(T)  is  defined  as  described.  Using  Proposition  2, 
the  entropy  condition  s(E)  implies  that  between  Z+  and  Z_  on  the  Hugoniot 
curve  with  center  Z+,  L(Z)  >  0  >  M(Z).  Property  i)  follows.  Property  ii) 
follows  from  i)  using  (2.7),  since  M  <  0  on  the  curve  (L  *  0).  The  "domain" 
condition  (1.13),  with  (1.11),  now  implies  l(x)  is  defined  as  described. 

The  existence  of  the  shock  layer  follows  from  i),  ii)  (see  Fig.  2). 

The  region  R  *  {( T , 0 ) |  x+  <  x  <  t  and  i(x)  <  0  <  m(x)  }  is  negatively 
invariant  under  the  flow  induced  by  (1.6).  Also,  any  trajectory  Z(£)  which 
starts  in  R  must  be  monotone  (  6  increasing,  r  decreasing)  and  tend  to  Z_  as 
S  ♦  — . 

Consider  a  vertical  line  segment  crossing  R, 

{(i,0)l  t  -  xQ  ,  Hr)  <  0  <  m(T)  )  ,  x+  <  xQ  <  x_  . 
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Figure  2.  Phase  portrait  of  (1.6)  for  a  compressive  shock 


A  point  on  this  segment  may  belong  to  one  of  three  disjoint  classes *  the 
forward  trajectory  of  (1.6)  starting  at  the  point  may 

a)  exit  the  region  R  on  the  curve  (L  «*  0) 

b)  exit  the  region  R  on  the  curve  (M  «  0) 

c)  not  exit  the  region  R  . 

In  the  last  case,  the  trajectory  must  tend  toward  Z+  by  monotonicity.  The 
first  two  classes  are  open  (continuity)  and  nonempty  (endpoints).  The  third 
class  is  therefore  also  nonempty  by  connectedness,  so  some  shock  layer  exists. 

He  turn  to  demonstrate  the  uniqueness  of  the  shock  layer.  First  we  claim 
that  any  trajectory  of  (1.6)  joining  Z+  and  Z_  must  lie  entirely  in  the  region 
R.  Indeed,  no  trajectory  can  approach  Z+  within  either  of  the  regions  (M  <  0 
and  I#  <  0)  or  (M  >  0  and  L  >  0).  Also,  the  region 

R,  *  {(t,0)  t  8.|  t  <  t  and  L>0  or  M>0} 

is  negatively  invariant.  So  our  claim  holds. 

The  characteristic  equation  for  (1.6)  at  the  critical  point  Z+  always  has 


one  positive  and  one  nonpositive  root: 


(2.9) 


0 


M^/y  -  k  MQ/y 

lt/x  I*e/X  -  tc 

-  k2  -  (MT/y  +  L6A>k  ♦  d'(T+)  -  m'(T+>)MeL9/yX 
The  constant  tern  is  nonpositive.  If  it  is  zero,  then  >  0  .  In  general, 
then,  Z+  is  a  saddle  point,  and  the  uniqueness  of  the  trajectory  approaching 
within  R  is  easily  established.  In  the  degenerate  case,  one  nay  construct  a 
one  dimensional  center-stable  manifold  in  a  neighborhood  of  Z+,  locally 
invariant  under  <1.6).  A  trajectory  starting  in  this  neighborhood  which  does 
not  lie  in  the  given  center-stable  manifold  must  eventually  (for  £  large) 
leave  any  sufficiently  small  neighborhood  of  Z+  (see  Kelley  [3]).  Therefore, 
any  two  trajectories  approaching  Z+  within  R  oust  lie  on  the  same  curve. 

We  have  established  the  existence  and  uniqueness  of  the  shock  layer  when 
T  >  t  and  the  entropy  condition  s(E)  is  satisfied.  If  T  >  t+  but  the 
condition  s(E)  fails  to  hold,  arguments  similar  to  those  above  shew  that  no 
trajectory  of  (1.6)  can  connect  Z_  to  Z+.  Theorem  2  (and  Theorem  1  for 
compressive  shocks)  follows. 

For  a  rarefaction  shock  Z+,  Z_  with  T_  <  t+  ,  s  <  0  ,  satisfying  the 
entropy  condition  s(E),  the  situation  is  as  follows:  L(Z)  <  0  <  M(Z)  for  Z 
between  Z+  and  Z_  on  the  Hugoniot  curve  with  center  Z_.  Hence  m(T)  <  *<t) 
for  t  <  x  <  t+  .  If  l'(t)  <  0,  then  the  region 

R  ■  {<t,9)|t_<t<t+  and  m(T)  <  0  <  !(t)  } 

is  negatively  invariant  and  the  existence  and  uniqueness  of  a  trajectory 
connecting  Z_  to  Z+  is  established  as  in  the  compressive  case  above.  However, 
it  no  longer  folows  from  (2.7)  that  l'(T)  <  0  if  t(T)  >  m(T)  .  But  we  may 
calculate,  using  (2.2), 

_“Lt<t,0)  -  eT<t,0)  +  p(x,0)  +  sM(t,0) 

■  eT(T,0)  +  p( T,m( t ) ) 

So  the  condition  eT<T,0)  >  0  ensures  that  t'(T)  <  0  for  all  t,  and 


Theorem  1  follows 


We  conclude  this  section  by  remarking  that  it  is  easy  to  use  Propositon 
2,  (2.7)  and  (2.8)  to  show  that  a  simple  jung>  solution  admits  a  purely  viscous 
shock  layer  (X  -  0)  if  and  only  if  the  entropy  condition  s(E)  is  satisfied. 
This  recovers  the  result  of  Liu  [4],  under  the  assumptions  (1.8)-(1.10),  which 
are  rather  different  them  those  Liu  imposed  (see  [S]). 


3.  A  shock  without  a  shock  layer 


In  thia  section  we  construct  explicitly  an  equation  of  state  satisfying 
assumptions  (1.8)— (1.10),  and  (1.14),  (1.15),  defined  in  a  suitable  region  of 
the  ( T, 8)  plane,  but  violating  the  restriction  6^(1, 6)  >  0,  such  that  a  (large 
amplitude)  rarefaction  shock  exists,  satisfying  the  entropy  condition  s(B), 
for  which  no  shock  layer  exists  if  \/v  is  sufficiently  large.  The  shock  layer 
does  exist  if  \/v  is  small,  and  we  stress  that  entropy  does  increase  across 
this  shock,  i.e.,  S(Z_)  <  S(Z+)  with  s  <  0  . 

Our  object  is  to  construct  Fig.  3,  representing  the  phase  portrait  for 

(1.6) .  The  idea  is  that  if  \/\i  is  large,  the  vector  field  (tH/u,  tL/X  )  for 

(1.6)  is  nearly  horizontal,  and  a  trajectory  leaving  Zg  must  hit  the  "hump"  in 
the  curve  (M  *  0).  Then  it  is  easy  to  show  that  no  trajectory  connects  z_  to 
Z+.  But  using  (2.6)  and  (2.7),  one  may  easily  check  that  the  entropy 
condition  s(E)  is  satisfied  for  the  jump  determined  by  Z+,  Z_,  and  s  <  0. 


Our  procedure  is  to  specify  selected  date  (in  particular  the  curve 
0  -  m(T)  on  which  M  •  0),  using  theraodynasdc  identities  to  consistently 
specify  the  equation  of  state*  Fix  Z0  ■  ( T0 , 00 )  in  the  first  quadrant.  We 
shall  specify  the  curve  (M  -  0)  as  0  -  m(T),  where  m(T)  satisfies  conditions 
(m1)~(m3)  below  (so  its  graph  appears  as  in  Fig.  3).  Later,  we  will  specify 
pQ  -  p(Tg,0g)  and  s.  Then  we  will  know  p  along  the  curve  (M  ■  0): 

(3. 1 >  pm<T>  5  p(T,n(T>)  -  pQ  -  a2(r  -  tq) 

Far  all  (t,0)  ,  we  specify  that 

(3.2)  0p0(r,0>  -  eT<T,0)  ♦  p(T,0>  • 

so  that  p.  -  1/0 2  ,  and  (1.10)  holds.  The  equation  of  state  for  p  is  now 
determined: 

(3.3)  p( T , 0 )  -  p  (T)  +  (i/m(T)  -  1/0) 

m 

Our  construction  will  be  complete  once  e0(T,0)  is  determined.  For  then, 
since  pQ,  a,  and  m(T)  determine  p,  eT(r,0)  is  known  from  (3.2).  Then  e^ 
and  e0  determine  e(T,0),  and  S(t,0)  is  determined  from  the  Gibbs  relation 
(1.2).  ((3.2)  ensures  that  the  expression  obtained  for  dS  can  be  integrated.) 

We  must  have  (e„)  -  (e  )_  *  -2/02  ,  so  e_  has  the  form 

OTTO  O 

e0(T,0)  -  e0<T0,0)  -  2(t  -  tq)/02 

We  specify  e„(T  ,0)  so  that  (1.8)  holds  in  the  region  of  interest.  It 
0  0 

suffices  to  define 

We>  “  8°TO(0"2  +  002) 

so  that  e.(x,0)  >  0  for  0  >  0,  0  <  T  <  40t  .  An  explicit  formula  for 

0  0 

e(T,0)  is 

17  T 

e(T,0)  -  80VV  ~  0  >e  _  /  V®( t)  dT 

T0 

+  (T  -  TQ) (2/0  -  P0  +  V*2(T  -  I0)>  +  e0 
We  must  now  show  that  m(T),  p0  and  s  may  be  specified  so  that  Fig.  3  is 
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(3.4) 


valid  and  (1.9)  holds,  along  with  (1.11),  (1.14),  (1.15).  Observe  that 

-sM(t,8)  is  independent  of  pg  and  ss 

-sM(x,0)  -  p(x,0)  -  p  (t)  -  1/m(x)  -  1/0 

m 

Using  (2.7),  we  must  have 

-L  ( T , 0 )  -  0pa  -  (-s)M  -  2/0  -  1/ra(T) 

-S  T  0 

Requiring  L(xo,0o)  ■  0  ,  the  quantity  <1/-s)L(x,0)  is  independent  of  pQ  and  s 
(recall  (1/-s)L0  -  eQ  ) s 

(3.5)  — L(T,0)  -  80Tn ( 0~2  -  0~2 ) 0  +  2( T  -  Tn)  -  /  1/m(T)  dt 

— S  0  0  u 


We  will  use  this  formula  to  manipulate  the  curve  (L  ■  0)  by  choice  of 
m(T)  .  Set  =  (k+1)xQ  .  We  shall  show  that  m(t)  can  be  chosen  so  that: 


i) 


L( t ,m( t ) )  >  0  for  x  -  xQ  >  0  small 


(3.6) 


ii)  L(t,0q+€1)  <  0  for  t3  <  T  <  T5  for  Bome  e1  >  0  • 
while  0Q  <  m(r)  <  0Q+e1  for  T4  <  T  <  T5 


iii)  L( t ,m( x ) )  >  0  for  T  large 

We  ask  that  m(x)  be  defined  and  smooth  for  <  t  <  4Ox0  ,  and  satisfy 

conditions  (m1)-(m3)  below: 

(ml)  m(xQ)  *  0o  and  m*(Tg)  *  ° 

Since  ( 1/-s)Lt(to,0o)  -  1/0Q  >  0  ,  <3.6i)  follows. 


(m2) 


m'(x)  <  0 
m(T) 
m'(T)  >  0 


for  tq  <  t  <  x, 


9g/6  for  x.|  <  T  <  T3 


for  x„  <  x  <  x.  ,  and 
3  4 


max  m( T )  *  0,+e.  ,  e.  >  0  small,  with  m(T_)  *  8_  *  m(T_) 
.  .  011  405 

VT<TS 


To  verify  (3.6ii),  we  estimate,  for  x^  <  x  <  x^  , 


-sL<T'W  <  CE1  +  2(T5  '  T0)/60  ‘  (T3  *  V6/0O 


Ce1  -  2Xg/0g  <0  if  e1  is  small. 
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(m3)  m(x)  -  0Q  for  x  >  xg 

To  verify  (3.6iii),  we  estimate,  for  T  >  x^  , 

_iut,e0>  >  2(,  -  T0)/«0  -  «vV6/eo  ♦  (T  -  t4)/o0> 

-  <(T  -  T4>  -  4(T4  -  T0)>/60  >0  if  T  >  T20  . 

Now:  Since  ( 1/-s)L0(x,8)  >  c  >  0  for  V/q  <  t  <  x4Q  ,  the  set  (L  ■  0) 

is  a  curve  0  *  4(  x )  which  is  bounded  above.  We  claim  that  £(x)  is  defined  for 

T.  <  x  <  t  .  Indeed,  L(t,0  )  <  0  for  t,  <  t  <  t_  ,  and  since  L  (x,0  /6)  >  0 
u  D  0  3  5  TO 

for  T  <  t3  (use  (2.7)  with  M  4  0),  we  have  L(t,80/6)  <  0  for  x  <  x^  , 
so  £(x)  is  bounded  below.  In  fact,  we  have  shown  that  £(x)  >  m(x) 
for  <  T  <  Tj.  . 

We  now  set 

t+  *  min{  T  |  t  >  ij  and  £(x)  *  m(x)  } 
t_  *  raax{  t  |  t  4  t  and  i(x)  *  m(x)  } 

Then  x  <  x  <  x  ,  and  since  m'(xft)  *  0  >  £'(xA),  xA  <  x  <  x,  .  We  set 

3  +  ZU  0  00“l 

0+  *  m(x+),  8_  ■>  m(x_),  fixing  2+  and  Z_,  and  completing  the  construction  of 
Fig.  3. 

It  remains  to  choose  pQ  and  s  so  that  p(x,0)  >  0  and  (1.9),  (1.14), 

(1.15)  hold  in  a  suitable  domain.  But 

<3*7>  Pt<t,S>  -  pT(x,0)  -  0pQ(x,0)2/e0(x,8> 

(One  calculates  p(x,S(x,0))  and  0  (x,S)  -  -S  /S «(x,0)  and  uses  (2.1), (2 .2).) 

T  T  TO 

So  (1.15)  implies  (1.9).  But  from  (3.2), 

(3.8)  Pt<t,0)  -  -s2  +  ( 1/m<  x) ) ' 

So  we  can  ensure  that  pt(t,0)  4  0  for  0  <  x  <  x4Q  if  we  choose  s2  sufficiently 
large. 

As  we  choose  pQ,  we  seek  to  fix  the  domain  of  states  in  the  (x,p) 
plane  so  that  flp  is  convex,  p  >  0  in  (1.14)  holds,  and,  pulled  back  to  the 
(x,8)  plane,  the  domain  O0  contains  (say)  the  half  strip 
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{(x,8)|  0  <  T  <  t  ,  fl  >  0/10  }  which  contains  all  our  constructions.  It 
40  0 

suffices  to  choose  ft  of  the  form 

P 

ftp  -  {<t,p)|  0  <  -t  <  x40  and  p  -  pQ  +  s2(t  -  xQ)  >  -1O/0Q  > 

with 

(3.9)  pfl  -  40tqs2  +  10/eQ 

To  check  (1.14),  compute 

e  (x,p)  -  e  (x,0)  +  eft(x,8)(-p  (T,0)/p  (x,0)) 

T  T  0  T  D 

2 

Using  (3.2)  and  (3.8)  one  may  verify  that  eT(x,p)  >  0  in  ftp  if  s  is 
sufficiently  large. 

Now  the  equation  of  state  is  completely  specified,  and  (1.8)— (1.15)  hold 
except  for  (1.12).  We  have  already  noted  that  the  entropy  condition  s(E) 
holds,  since  i(x)  >  m(x)  for  t_  <  t  <  x+  . 

The  penultimate  stage  of  our  analysis  is  to  show  that  for  X/ji 
sufficiently  large,  there  is  a  trajectory  of  (1.6)  leaving  ZQ  and  intersecting 
the  curve  (M  »  0)  before  the  "hump".  Since  m'(xQ)  >  i'(xQ),  ZQ  is  a  saddle 
point  for  (1.6)  (see  (2.9)).  If  k  is  the  positive  characteristic  value  for 
(1.6)  at  Zq,  the  the  characteristic  vector  (x  -  XL^/X , xL^/X )  may  be  shown  to 
have  positive  components.  Let  Z(£)  be  the  solution  leaving  Zq  along  this 
vector.  So  long  as  Z(5)  remains  above  the  curve  (M  =  0),  x  is  increasing  on 

Z(C)  and  we  may  write  0  *  z(x)  along  Z(£).  Then 

dz  _  mL(x,z(x) ) 
dx  "  XM(x,z(x)) 

It  is  clear  that  z(t)  is  defined  and  increasing  for  xQ  <  x  <  t#,  where 
x#  ■  min{  x  |  x  >  xQ  and  L(x,0Q)  ■  0  ) 

Note  that  x#  <  x3«  We  will  "channel"  the  curve  z(x)  under  the  curve  (L  -  0) 
by  using  a  thin  box  with  a  corner  cut  outs  Define 

R*  ‘  (VT3)X<VE2'W  ^  <T0'T*)X<V£2'e0^1> 

Here  >  0  is  chosen  so  that  R*  is  contained  in  the  region  (M  >  0).  Now,  the 


trajectory  Z(£)  must  enter  R*  on  the  cut-out  pert  of  the  boundary  with 
9  >  9Q,  i.e.,  it  nut  enter  at  a  point  (T^-tyf ^ )  for  tq  <  T  <  t#,  or  at  a 
point  (t#,8)  for  8^  <  8  < 

Now  sup  |L/M)  is  finite.  If  X/u  is  sufficiently  large  in  R* ,  we  nay 
R* 

obtain 

I  H  I  <  “in {  1/( t3”t0 > •  e2/(t3-t#>} 

so  long  as  Z(£)  is  in  R*.  This  inplies  that  the  trajectory  Z(  t)  oust  leave  R* 
on  its  right  boundary,  }x  (®g"e2  '®0+et^ '  wh^ci*  *■*  contained  in  the  region 
(I*  <  0  and  M  >  0).  Past  this  point,  8  decreases  and  t  increases  along  Z(  £) , 
so  since  m(T)  acheives  the  value  8^+6^,  Z(  5)  oust  intersect  the  curve  (M  »  0) 
at  a  point  (t,8),  where  >  t  >  >  t 

We  nay  now  easily  construct  a  positively  invariant  region  containing  Z_ 
and  disjoint  fron  Z+,  so  that  no  trajectory  can  connect  Z_  to  Z+.  The  region 
is  bounded  on  the  right  by  the  line  t  •  T  for  8  <  8  (M  <  0  here),  from  above 
by  the  curve  8  *  z(t)  for  tq  <  T  <  T  (trajectories  cannot  cross),  and  by  the 
curve  9  *  min{l(T),m(T>}  for  T  <  Tg  (the  vector  field  either  points  down,  or 
to  the  left  at  a  point  where  t'(t)  <  0,  as  in  §2).  Our  example  is  complete. 
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